Effect of disorder on the electronic properties of graphene: A theoretical approach 
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In order to manipulate the properties of graphene, its very important to understand the electronic 
structure in presence of disorder. We investigate, within a tight-binding description, the effects of 
disorder in the on-site (diagonal disorder) term in the Hamiltonian as well as in the hopping integral 
(off-diagonal disorder) on the electronic dispersion and density of states by the augmented space 
recursion method. Extrinsic off-diagonal disorder is shown to have dramatic effects on the two- 
dimensional Dirac-cone, including asymmetries in the band structures as well as the presence of 
discontinuous bands (because of resonances) in certain limits. Disorder-induced broadening, related 
to the scattering length (or life-time) of Bloch electrons, is modified significantly with increasing 
strength of disorder. We propose that our methodology is suitable for the study of the effects of 
disorder in other 2D materials, such as a boron nitride mono layer. 



PACS numbers: 73.22.Pr, 61.48.Gh 

I. INTRODUCTION 

Graphene, a two-dimensional allotrope of carbon, plays 
a central role in providing a basis for understanding the 
electronic properties of other carbon allotropes. Being 
one of the thinnest and the strongest material ever mea- 
sured, graphene has attracted the attention of the ma- 
terials research community 1 in the recent past. One of 
the most interesting aspects of graphene is that its low 
energy dispersion closely resembles the Dirac spectrum 
of massless fermions. This particular type of dispersion 
provides a bridge between condensed matter physics and 
quantum electrodynamics (QED) for massless fermions. 
Of course in graphene, the Dirac fermions move with a 
much smaller speed. 

Because of its unusual electronic and structural flexi- 
bility, properties of graphene can be controlled chemically 
or structurally in many different ways. For example, de- 
position of metal atoms 2 on top of the graphene sheet, 
incorporating other elements like boron and nitrogen 3 
randomly in the parent structure, cither interstitially or 
substitutionally and using different substrates. 4 Because 
disorder is unavoidable in any material, there has been 
an increasing interest in understanding how disorder af- 
fects the physics of electrons in graphene. 5 Disordered 
graphene based derivatives can probably be referred to 
as functionalized graphene suitable for specific applica- 
tions. " Graphene paper 1 ' 6 is a spectacular example of 
how important such functionalization could be. 

There can be many different sources of disorder in 
graphene including both intrinsic as well as extrinsic. 



Intrinsic sources may include surface ripples and topo- 
logical defects. Extrinsic disorder comes in the form of 
vacancies, adatoms, quenched substitutional atoms and 
extended defects such as edges and cracks. Another way 
of introducing disorder is by ion-irradiation that pro- 
duces complex defect structures in the graphene lattice. 7 
Graphene in an amorphous form may increase the metal- 
licity too. 8 

To have a theoretical description of graphene's elec- 
tronic structure, one may begin with the Kohn-Sham 
equation and a tight-binding representation whose ba- 
sis is labeled by the sites of the underlying Bravais lat- 
tice. Disorder may enter the matrix representation of 
the Hamiltonian in two ways : vacancies, dopants and 
adatoms predominantly cause a random change in the 
local single-site energy (disorder in the diagonal terms) 
but through the overlap such defects modify the hop- 
ping integrals between different sites (disorder in the off- 
diagonal terms) causing an effective random change in 
the distance or angle between the bonding orbitals. Thus 
diagonal and off-diagonal disorders simultaneously occur 
and are correlated. Model calculations which take them 
to be independent are qualitatively in error. As far as 
diagonal disorder is concerned, it acts as a simple chem- 
ical potential shift of the Dirac fermion i.e. shifts the 
Dirac point locally. Theoretical study of such disorder 
is rather simple and has indeed received attention and 
success, reported in literature. 5 ~ 9 A proper inclusion of 
off-diagonal disorder, on the other hand, is non-trivial 
and requires more sophisticated approaches. 

Till date there have been numerous attempts at study- 
ing the effects of disorder in graphene. 5-9 Among others, 
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FIG. 1: (Color online) (Left) The standard honeycomb lattice 
with a basis of two atoms per unit cell. (Right) The under- 
lying rhombic Bravais lattice which becomes the honeycomb 
lattice when a pair of atoms decorate each site. 



the methods used to study effects of disorder included 
the averaged t-matrix approximation (ATA) 9 and the co- 
herent potential approximation(CPA). 10 The first one is 
not self-consistent and hence inaccurate. The latter is 
a single-site mean field approximation with all its atten- 
dant problems. 10 Several others have used exact diago- 
nalization of huge clusters and the real-space recursion 
of Haydock et al. 9,11 Both these techniques actually cal- 
culate the density of states (DOS) for specific configura- 
tions of the system followed by direct averaging over a 
large number of configurations. Since each of the con- 
figurations has periodic boundary conditions, the aver- 
aged spectral function is always a collection of delta func- 
tions and the disorder induced life-time effects cannot be 
probed. The recursion on the lattice probes mainly the 
real-space effects of disorder. 

From the theoretical perspective, dealing with disor- 
der has had a long history. As mentioned earlier, one 
of the most successful and frequently used approaches is 
the single-site, mean field CPA. 10 However, as the name 
itself suggests, it is a single site approximation and can- 
not adequately take into account the effects of correlated 
configuration fluctuations. In particular the CPA is in- 




FIG. 2: (Color Online) Nearest neighbor overlaps on the 
rhombic lattice 



accurate at low dimensions. In one dimension it is shown 
to be inadequate by Dean 12 some time ago. Among the 
hierarchy of the generalizations of CPA, only a few ap- 
proaches have maintained the necessary Herglotz ana- 
lytic properties and lattice translation symmetry of the 
configuration averaged Green's function. These include 
the non-local CPA, 13 the special quasi-random structures 
(SQS), 14 the locally self-consistent multiple scattering 
approach (LSMS), 15 and the three methods based on the 
augmented space formalism proposed by one of us 16 : 
the traveling cluster approximation (TCA), 17 the itiner- 
ant coherent potential approximation (ICPA), 18 and the 
augmented space recursion (ASR). 19 Over the years ASR 
has proved to be one of the most powerful techniques, 
which can accurately take into account the effects of cor- 
related fluctuations arising out of the disorder in the lo- 
cal environment. This is reflected in a series of studies 
in the past e.g. the effects of local lattice distortion as 
in CuBe, 20 short-range ordering due to local chemistry, 21 
the phonon problem 22 with essential off-diagonal disor- 
der in the dynamical matrices, and electrical and thermal 
transport properties 23 in disordered alloys. 

In this communication, we present a theoretical tight- 
binding model to study the effects of disorder in 
graphene. Disorders studied were mainly of two forms 
: substitutional disorder 24 ' 25-26 and vacancies. 9,27 Un- 
like earlier models, both the diagonal and off-diagonal 
disorders are included on the same footing. The present 
formalism is based on the augmented space recursion. 19 
Although recursion has been used to study graphene be- 
fore, we want to emphasize that in all those applications 
recursion was carried out on a Hilbert space H spanned 
by the tight-binding basis representing the Hamiltonian. 
In augmented space recursion, we recurse in the space of 
all possible configurations which the Hamiltonian may as- 
sume in the disordered system. For a homogeneously dis- 
ordered binary alloy, this configuration space is isomor- 
phic to that of a spin-half Ising model. The augmented 
space theorem 16 then connects configuration averages to 
a specific matrix element in that space of configurations. 

The novel approach in this work is that we shall make 
use of the translation symmetries in augmented space (for 
homogeneous disorder) to carry out recursion in recipro- 
cal space. This will directly give us the spectral function 
from which we extract the 'fuzzy' band structure. The 
inclusion of the effects of configuration fluctuations of the 
immediate environment gives us self-energies which are 
strongly k dependent, unlike the CPA. In order to make 
a systematic study, we present results for combinations 
of both strong and weak diagonal and off-diagonal disor- 
der. The combined effects show dramatic changes in the 
location and topology of the Dirac-like dispersion and the 
DOS. Special emphasis has been given to the non-trivial 
inclusion of off-diagonal disorder, in which case the aver- 
aged Bloch spectral function comes out to be significantly 
broadened, multiply peaked, and asymmetric in certain 
limits where the presence of resonances leads to discon- 



3 



tinuous dispersion. The interesting interplay of the two 
kinds of disorder on full-widths at half maxima (FWHM) 
(related to the life-time of Bloch electrons in a disordered 
system) is also shown. 

The rest of the paper is organized as follows. In Sec. II, 
we introduce the basic formalism. Sec. Ill is devoted to 
results and discussions. Concluding remarks are present 
in Sec. IV. 



For a system with substitutional disorder, the most 
general statement we can make is that the occupation 
of the lattice sites in each inequivalent sublattice can be 
different. For binary disorder in both the sublattices, we 
may introduce two random occupation variables n R and 
n R associated with the sublattices / and II such that, 



1 if R e A with probability 
if R G B with probability x B 



and 



II. FORMALISM 



The most general tight-binding Hamiltonian for elec- 
trons in graphene can be represented as, 

^ = EE {^SrwSss^+V^'T^'}, (1) 

Ra a R'a s , 

where R, R' denotes the position of the unit cell of the 
lattice, a s denotes the a-th atom on the s-th sublattice. 
The actual atomic position is R + ( as , where ( Us is the 
position of the a-th atom on the s-th sublattice. e^ s 
is the on-site energy describing the scattering properties 
of the atomic potential at R + ( as , and V RR J 3 ' is the 
hopping integral between R + ( a and R' + ( a . P and T 
are the projection and transfer operators in the Hilbcrt 
space spanned by the tight-binding basis \Ra s ). 

The above Hamiltonian H describes electrons in the 
original honeycomb lattice of ion-cores, as shown in the 
left panel of Fig. 1. The two inequivalent sublattices 
(shown by red and blue spheres) are distinguished from 
each other. The underlying Bravais lattice is the rhombic 
lattice shown in the right panel of Fig. 1. Looking at 
Fig. 2 we can simplify Eqn.(l) further and write the 
Hamiltonian elements as 2 x 2 matrices : 



H = Y,iR p *+ E ( 2 ) 

R R^R' 

where e R and ¥_ , > instead of being scalar for a single- 
band problem, are now 2x2 matrices given by, 



= ei t /f 

~ R \ t £ 2 I =01 ^02 I t t' 



^03=^04=1 I )' ' 



where e\ and £2 are the on-site energy on the two sub- 
lattices, t and t' are the nearest neighbor and the next 
nearest neighbor hopping energies. V_ QI are the hopping 
matrices between the central site and its four neigh- 
boring sites / (in the rhombic lattice) as shown in the 
right panel of Fig. 1. Because the next nearest neighbor 
hopping t' is usually very small compared to t, we shall 
treat the disorder effects only in the nearest neighbors. 



n R 



1 if R e C with probability xc 
if R e D with probability x B 



where A, B are the two types of atoms randomly occupy- 
ing sublattice / and C, D are those occupying sublattice 
II. 

The diagonal term £ R for such a binary distribution 
can be written as, 



e A tAC \ 1 11 

JI n R n R 



t AC e J c J 



+ 



= E 1 +E 2 n I R + E 3 71% + Ei n^njl, 
where 



(4) 



Ex = 
E 3 — 



£3 ^BD 
tBD £ B 

t& 

5£ 2 



Eo = 



E A = 



tw 



t< 3 ) 



(5) 



with Sei — £ A — £ B ; with fc 2 = £q — £ B : = tAD — 
tBD, i (2) = tBC-tBD and = t AC + t BD -t AD -t BC - 

Similarly the off-diagonal term V RR i in Eq. 3 can be 
expressed as (assuming t' = 0), 

Z01 = Z02 = V 1 + V 2 n R + V 3 n R \ + V 4 n£n# , (6) 
where 



V* 





t BD 







v 2 



v 4 





t« 

f <°3, Si- (A 



V_ Q3 (= Y_ m ) are just the transpose of the above matrix 
V 01 . Various t a ps in the above sets of equation are the 
hopping energies between various atom types (a = A, B 
and /? = C,D) at two sublattices / and II respectively. 

Next we proceed to calculate the configuration aver- 
aged Green function ( or the Bloch spectral function) 
in reciprocal space. We shall generalize the augmented 
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space formalism (ASF) developed earlier in reciprocal 
space. 28 The ASF has been described in great detail 
earlier. 29 We shall indicate the main operational results 
here and refer the reader to the above monograph for 
further details. The first step is to associate with n R and 
n R two operators N R and N R such that their spectral 
density is the probability density of the random variables. 
For binary random variables, we have : 



R \ yJx A X B X A J 

Finally, according to augmented space theorem, 16 the 
configuration average of any function of {n R ,n R } can be 
written as the matrix element, in configuration space, of 
an operator which is the same functional of {N R7 N R }. 
The augmented space Hamiltonian is built up from Eqns. 
(4) and (6). 



< G(k,z) >= ({0}®k|(zI-H)- 1 | {0}), (9) 

| k ® {0}) is an augmented space state in the reciprocal 
space given by, 



|k®{0}> = ^£y ik - R |i2®{0}>, (10) 

V JV R 

| i?®{0}) is an enlarged basis which is a direct product of 
the Hilbert space basis {R} and the configuration space 
basis {4>r}. The configuration space $ = FJl <j) R , takes 
care of the statistical average, is of rank 2 M for a system 
of M-lattice sites with binary distribution. 

The recursion follows as a three step generation of a 
new basis {\n >} : 



H = {EiT+ E 2 N r + EzNH + E A N R ® Ntf} <8> P R 

R 

+ ]T ]T { V ^+ V ^R + V ^R' + V ^R ®N I R 1 ,}® T RR , 

R R> 

with 

= *a VV + Z/3 PR ^ + V^( T R U + T R M )' (8) 

(X = I or II) 

The configuration averaged Green's function in the re- 
ciprocal space is thus a matrix element of an augmented 
resolvent given by, 



|l) = |k®{0}> |0>=0 

\n + l)=H\n)-a n \n)-l3 2 n _ l \n-l) 

j 2 n , (n|n) 
a " (k) = W ^ ^ (k) = (n-l|n-l) 

The ASR gives the configuration averaged spectral 
function as a continued fraction : 



< G(k,z) > = 



z - ai(k) 



/?i 2 (k) 



z-E Q {k)-Y,{z,k) 



(11) 



z - a 2 (k) 



T(z, k) is a continued fraction terminator as proposed 
by Beer and Pettifor. 30 The spectral function peaks are 
decided by 5Re{E(i?,k)} and the imaginary part of E 
gives the width related to the disorder induced lifetimes. 

The configuration averaged Bloch spectral function is 
given by, 

«A(k, £)>=-- lim 3?m{< G(k, E + iS) >} (12) 

7T (5->0+ 



The configuration averaged density of states (DOS) is, 
< n{E) >= — — / dk < A(k, E) > (13) 

"BZ JBZ 

The electronic dispersion curves are obtained by nu- 
merically calculating the peak E-position of the spectral 
function. The full-widths at half maxima (FWHM) are 
also calculated from the disorder broadened Bloch spec- 
tral function. 
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FIG. 3: Local density of states for pure graphene and graphene with a single impurity. The top panel in the middle is the DOS 
for pure graphene. Left and right panels show the local DOS with single impurity only on sublattice I or II and both I &: II 
respectively. The panel from top to bottom are the results with increasing strength of the impurity potential <5E=£i mp -£host 



III. RESULTS AND DISCUSSION 



A. Impurities in Graphene 



In the following subsections, we shall present our re- 
sults for graphene with impurities, vacancies, diagonal 
disorder alone, and with the simultaneous presence of 
diagonal and off-diagonal disorder. The effects of var- 
ious strengths of impurity potentials on two inequiva- 
lent sublattices will be shown via changes in the shape 
of the DOS. The changes in the topology of Dirac-cone 
dispersion, disorder-induced FWHM and the DOS will 
be shown for various strengths of diagonal disorder. In 
the most general case of diagonal and off-diagonal dis- 
order, we consider three interesting limiting cases: (i) 
strong diagonal and weak off-diagonal disorder (ii) strong 
off-diagonal and weak diagonal disorder and (iii) strong 
diagonal as well as off-diagonal disorder. The interest- 
ing interplay between these different kinds of disorder in 
graphene reveals a discontinuous type of band near the 
T-point in the third limiting case. 



In Fig. 3, we display the DOS with different strengths 
of the single impurity potentials on different inequivalent 
sublattices. The top figure in the middle panel is the 
DOS for pure graphene. The left and right panels show 
local DOS with a single impurity put on the sublattice I 
or II and both respectively. The strength of the impurity 
potential (relative to the host lattice) increases from top 
to bottom panels (i.e. 5F, = e imp - £ host — 0.4, 0.7 and 
1.0). All these calculations are done with a fixed hop- 
ping parameter t=l. We notice changes in the shape of 
the hump and the van-Hove singularities as the strength 
of the impurity potential increases. Although the effects 
are small, but are clearly visible for the case of <5E=1.0, 
where the local environment around the impurity site 
feels the strongest scattering. With the introduction of 
the impurity, the symmetry of the DOS around the Dirac 
point is lost. At these impurity levels, both the left and 
right panels show the formation of an impurity peak near 
the upper band edge. With increasing disorder this im- 
purity peak moves into the band and disappears. Again, 
at these strengths there is no perceptible changes to the 
linear structure of the Dirac point. Similar results have 
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Energy (E) 

FIG. 4: Local density of states for graphene with a single 
vacancy. The vacancy site is modeled by a site with a large 
repulsive local potential. Technically we take <5E = £i mp — 
£host- The figures show consecutive situations with increasing 
SE as we go from top to bottom. 

been obtained previously for such models of impurities. 
This provides the correctness of our new formulation. 



B. Single vacancies in graphene 

An extension of the single impurity is the single va- 
cancy in the model. The vacancy is modeled by a site 
with a large repulsive local potential. This is shown in 
Fig. 4. Notice that, with increasing SE, the rightmost 
impurity peak at around the top band-edge moves into 
the band. Most of the changes occur around the Dirac 
point at E=0. At SE = 2 > t(= 1) the symmetry of the 
Dirac point gets broken with the appearance of another 
mild peak. The origin of this 'zero mode' peak has been 
extensively discussed by Pereira et al. 31 in detail. This 
peak grows with increasing disorder until at SE ~ 30, 



where it takes the form of a delta function like structure 
at the Dirac point. The bottom panel shows the result 
for an ideal vacancy, where Uj = for t connecting the 
vacancy to the graphene lattice, or in other words a com- 
pletely inaccessible 'hard' vacancy. The symmetry of the 
Dirac point is restored and the vacancy peak sits exactly 
at the Dirac point. This is exactly the same behavior 
reported by Pereira et al. 31 who used either the CPA or 
direct real space recursion. 



C. Diagonal disorder 

First we shall take up purely diagonal disorder prob- 
lems : those problems which can be taken up by ear- 
lier suggested methodologies. Of course, our augmented 
space recursion in reciprocal space give us additional in- 
formation about the disorder induced life-times of the 
Bloch states. In Fig. 5, we display the configuration 
averaged Bloch spectral function (upper panels) (given 
by Eq. (12)) and the corresponding dispersion Energy 
vs. k (lower panels) along r — K — T symmetry line for 
three different alloys A XA B XB (xb = 5%, 15% and 25%) 
with various diagonal disorder strengths. The two sublat- 
tices I and II are homogeneously disordered, such that 
xa = xc, xb = xd and e\ — Sq , e B = ejj . For each 
alloy case, the panels from top to bottom indicate the 
results with increasing strength of diagonal disorder (i.e. 
5E= e A - e B = 0.1,0.4,0.7 and 1.0). The hopping inte- 
gral t is chosen to be 1 here, so there is no off-diagonal 
disorder. For each disorder strength SE in a particular 
alloy, the upper panels show the averaged Bloch spec- 
tral function at five k-points along the high symmetry 
direction T— K. The corresponding Dirac dispersions are 
shown in the lower set of panels along r — K — T line. 
The first thing to note is that the spectral function mod- 
ifies quickly from sharp near <5-functions to Lorentzian 
shapes with increasing disorder strength SE as well as 
increasing alloy concentration xb- In addition, the func- 
tion gets more and more asymmetric with increasing SE. 
Such asymmetries can be described as a tendency of more 
scattering to occur near the resonance energies around T. 
In other words, line shapes around T tend to have a weak 
second peak or wide tail over the resonance region. For 
the present diagonal disordered case, the Dirac point is 
simply shifted by an average energy (e) = xa^\ + xb^J 

The corresponding total DOS for the same set of dis- 
order strengths (SE) and the alloy concentrations (x) are 
shown in Fig. 6. The individual projected DOS on the 
two sublattices I and II in this case are same, because 
we have maintained uniform diagonal disorder on both 
the sublattices. However, the present theory is equally 
capable of treating the two sublattices differently with a 
different nature of disorder on them. In that case, the 
two inequivalent sublattices will have different projected 
quantities. Looking at Fig. 6, one can notice an exactly 
similar shift of the Dirac point (to the average (e)) in 
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FIG. 5: (Color online) The averaged spectral functions (upper set of panels) and the complex bands (lower set of panels) near 
the Dirac point. These are all for pure diagonal disorder at three different alloy compositions (left to right) and four disorder 
strengths <5E (top to bottom). The (red) error bars show how the disorder induced lifetimes vary across the samples. 
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FIG. 6: Total DOS for the same set of disorder strengths <5E for three alloys K XA Q XB as in Fig. 5. Due to homogeneous 
diagonal disorder on both the sublattices I and 77, the individual projected DOS on them are same in this case. 



the DOS as shown in the dispersion. The disorder effects 
are pronounced around the Dirac-point energy (e) and 
get milder around the hump below SE = 0.7. Above this 
disorder strength, the left band edge starts to show up 
extra features with a dip at around E = — 2, (as shown in 
the bottom panel for all the three alloy concentration). 
The results are qualitatively similar to the CPA works 
done earlier 31 but differ in quantitative details. 



t AC = 1.0, t BD = 0.9 and t AD = t BC = 0.95. 

• Weak diagonal and strong off-diagonal disorder; 
with parameters SE = 0.1, t A c = E0,tBD = 0.5 
and t A D =tsc = 0.75. 

• Strong diagonal as well as strong off-diagonal dis- 
order; with parameters SE = 1.0, t A c = 1-0, tBD = 
0.5 and t AD = t BC = 0.75. 



D. Off-diagonal disorder 

We now turn to the cases with off-diagonal disorder. 
Such problems cannot be dealt with within the CPA. 
Also, direct calculation of the averaged spectral functions 
and disorder induced lifetimes is also not feasible with 
other techniques and the strength of the ASR comes to 
the fore. In addition we should note that in our model, 
diagonal and off-diagonal disorders are correlated : e.g. 
if the atom A occupies the site i with probability x A and 
atom B occupies the site j with probability Xb, then tij 
has to be t A s with probability 1. Although the present 
theory is equally capable of investigating other interest- 
ing cases (e.g. inhomogeneous disorder, pseudo-binary 
type disorder etc.), here we have chosen to explore three 
cases which should reflect the behavior of a variety of the 
realistic materials. The three cases are: 

• Strong diagonal and weak off-diagonal disorder; 
with parameters <5E = e A — e B = — e D ! = 1.0, 



The results for these three cases are shown in the top, 
middle and bottom panels of Fig. 7 respectively, for the 
same three alloys A XA B XB as before. Other details are 
same as in Fig. 5. Notice that unlike the diagonal dis- 
ordered case, effects of both diagonal and off-diagonal 
disorder are much more dramatic. In addition to highly 
asymmetric nature, the Bloch spectral function is found 
to have a double peaked structure in the extreme case 
of strong diagonal and off-diagonal disorder (shown in 
the bottom panels). Such doubly peaked line shape 
introduces extra discontinuous bands in the dispersion 
curve. Such a structure had been seen before in phonon 
problems 22 which also have intrinsic off-diagonal disor- 
der in the dynamical matrices. There it arose because of 
resonant modes. Here too we shall give a similar explana- 
tion. These dispersion at resonance have relatively large 
FWHM's and it will be interesting to choose a realistic 
material of similar disorder properties and investigate the 
experimental outcome. 

Figure 8 shows the sublattice projected DOS for the 
same three limiting cases for the three alloys as above. 
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The solid blue and the dashed red lines indicate the pro- 
jected DOS on the sublattices / and II respectively. Be- 
cause of the random hopping (off-diagonal) interaction 
in this case, the two sublattices acquire different envi- 



ronment around it, and hence possess different projected 
quantities on them. As expected, the DOS in these cases 
have large smearing. The effective environment around 
the two sublattices is maximally different from each other 
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FIG. 8: (Color online) Sublattice projected DOS for the same three set of diagonal+ off-diagonal disorder strengths for three 
alloys A XA B XB as in Fig. 7. The projected DOS on the two sublattices nj and nu, in this case, are different due to the obvious 
reason arising from different random tij interactions. 



in the extreme case of both strong diagonal+off-diagonal 
disorders (as shown in the bottom panels), as evident 
from the large difference between their projected DOS. 
Interestingly, for this particular case, the consequence of 
the discontinuous bands in the -ve energy range (bottom 
panels of Fig. 7) shows up via a dip in the DOS along 
with a much larger smearing. Apart from this extreme 
case, the Dirac point for all the other cases has moved in 
exact accordance with that of the band shift as in Fig. 7. 
The topology of the DOS on the two sides of the Dirac 
point are very different from each other specially in the 
case of strong diagonal and off-diagonal disorder (bottom 
panels) . In totality, the effects of off-diagonal disorder is 
very different from that of diagonal disorder (as a com- 
parison between Figs. 6 and 8 will show). Treatments of 
off-diagonal disorder is straight-forward and accurate in 
the ASR formalism. 



of off-diagonal disorder in describing the spectral proper- 
ties of graphene. We show how the topology of the Dirac 
dispersion, and the location of Dirac point change with 
the strength of disorder and impurity concentration. In- 
terestingly, the dispersion in case of strong diagonal and 
off-diagonal disorder tends to have an extra discontin- 
uous band which is rather uncommon in the graphene 
fermiology with simple disorder. As such we propose to 
verify such effects in the electronic dispersion by setting 
up an experiment on a similar realistic graphene system, 
where both the diagonal and the off-diagonal disorder 
are strong. We believe that such a study may provide a 
deeper insight into the physics and materials perspective 
of graphene. Finally, we want to state that our formu- 
lation is quite general and can be applied to the case of 
other 2D materials, e.g., BN in presence of disorder. 



IV. CONCLUSION 

We present a theoretical model to study the effects 
of diagonal and off-diagonal disorder in graphene on an 
equal footing. To our knowledge this is the first theoret- 
ical framework to reliably take into account the effects 
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